A random phase approximation study of one-dimensional fermions after a quantum 
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The effect of interactions on a system of fermions that are in a non-equilibrium steady state due 
to a quantum quench is studied employing the random- phase-approximation (RPA). As a result 
of the quench, the distribution function of the fermions is highly broadened. This gives rise to 
an enhanced particle-hole spectrum and over-damped collective modes for attractive interactions 
between fermions. On the other hand, for repulsive interactions, an undamped mode above the 
particle-hole continuum survives. The sensitivity of the result on the nature of the non-equihbrium 
steady state is explored by also considering a quench that produces a current carrying steady-state. 
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I. INTRODUCTION 

Recent remarkable experiments^ with cold atoms have 
motivated an explosion of theoretical interest in the area 
of non-equilibrium quantum dynamics with a focus on 
addressing fundamental questions about thermalization, 
chaos and integrability, issues that are very relevant to 
these experimental systems.^ Without many general re- 
sults on generic non-equilibrium phenomena, the analysis 
of specific, tractable models is a common way to make 
progress. One hopes clues gathered from these specific 
systems will lead to more general predictions. 

One-dimensional (ID) systems are where much of the 
theoretical work has taken place since a wide array of 
tools is available for investigating dynamics. An inter- 
esting class of these systems is integrable models, where 
conserved quantities tightly constrain the time evolu- 
tion. While a consensus lacks on a rigorous definition 
of quantum integrability^ progress has been made us- 
ing many quantum models satisfying classical notions of 
integrability. Fruitful studies have investigated dynam- 
ics of Bethe-ansatz solvable modelsj^ but the simplest 
integrable models are the quadratic ones. These effec- 
tively non-interacting theories, including those consid- 
ered in this paper, allow for exact analytical treatment 
of the non-trivial dynamics. In ID, efficient numerical 
studies are also possible with the time-dependent density 
matrix renormalization group (tDMRG),^^'^^ and exact 
diagonalization studies of finite systems. 

Some of the analytical and numerical studies have 
revealed that ID systems after a quantum quench of- 
ten reach athermal steady states which can be char- 
acterized by a generalized Gibbs ensemble (GGE) con- 
structed from identifying the conserved quantities of the 
system i^iiiiSiii^ii^ There are also many counter-examples 
where such a description fails, as not all physical quanti- 
ties can be described using the GGE<^iiii2r— 

One important question concerns the stability of these 
athermal steady states generated after a quantum quench 
to other perturbations such as non-trivial interactions 



that introduce mode-coupling and/or the breaking of 
integrability. Precisely this question was addressed re- 
cently in Ref [2ll In particular an initial interaction 
quench in a Luttinger liquid gives rise to an ather- 
mal steady state characterized by new power-law expo- 
nentSf^iiiSi which can also be captured using a GGE. 
The effect of mode-coupling arising due to a periodic 
potential on this non-equilibrium state was studied in 
Ref. I21I using perturbative renormalization group. The 
analysis revealed that infinitesimally small perturbations 
can generate not only an effective-temperature but also 
a dissipation or a finite lifetime of the bosonic modes. 
While the appearance of an effective-temperature, al- 
though highly non-trivial in itself, can be rationalized on 
the grounds that a system after a quench is in a highly 
excited state, and that interactions between particles will 
somehow cause the system to "thermalize" , the appear- 
ance of dissipation is an unexpected and non-trivial re- 
sult. Thus one of the motivations of the current paper 
is to identify other physical situations where this dissi- 
pation might appear, and to try to investigate the phys- 
ical mechanisms that could be behind it. Due to the 
close parallels between interacting bosons and fermions 
in ID, a natural candidate for analyzing this question is a 
one-dimensional system of free fermions that is in a non- 
equilibrium steady state after a quench. We analyze the 
effect of weak interactions on this system by employing 
the random-phase-approximation (RPA). 

In equilibrium, ID systems are the ideal playground 
for invoking the RPA. It is an exact low-energy treat- 
ment of weak interactions in IDi^^— In particular, by 



,26^ 



one re- 



applying RPA to a ID system of electron 
covers the standard bosonization of the model, described 
by a Luttinger liquid.^ Note that this direct equivalence 
only holds for the long wavelength properties, while other 
excitations require more sophisticated methods such as 
bosonization. While the accuracy of RPA in ID is known 
in equilibrium, its applicability out of equilibrium is not 
guaranteed. In the type of quench problem considered 
here, the initial state has nonzero overlap with excited 
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eigenstates of the Hamiltonian generating time evolution. 
It is far from certain that a low-energy description cap- 
tures all the important physics. While this caveat leads 
to intriguing, unanswered questions, we will use in this 
paper RPA as an approximation scheme and will not ad- 
dress the deeper question of its potential breakdown out 
of equilibrium. 

In this paper, we thus apply the RPA to a non- 
equilibrium state in the XXZ spin-chain. This state 
is prepared as follows. The system is initially in the 
ground state of an exactly solvable Hamiltonian Hi. We 
choose two different models for Hi, one corresponding to 
the transverse-field Ising model with the magnetic field 
tuned to the critical value where the spectrum is gapless, 
and the second is the same as above but with an addi- 
tional Dzyaloshinskii-Moriya interaction added. A quan- 
tum quench is then performed by switching off the field 
and changing the exchange anisotropy so that the time 
evolution is due to the XX model. Since this model is de- 
scribed by free fermions, at long times after the quench, 
the system reaches an athermal steady state character- 
ized by a GGE. For Hi that has Dzyaloshinskii-Moriya 
interactions, the steady state is qualitatively different in 
that it carries a net current. We then ask how these 
athermal steady states are affected by weak Ising inter- 
actions of the XXZ chain ( S^S^j^-^) which are as- 
sumed to have been switched on very slowly. The effects 
of the Ising interactions are treated using RPA. 

We demonstrate the existence of a single undamped 
collective mode for repulsive interactions ( Jz > 0) which 
is qualitatively similar to the predictions of the RPA in 
equilibrium, but with some quantitative changes to the 
mode-velocity. On the other hand for attractive interac- 
tions [Jz < 0), no undamped modes are found for both 
the athermal states that have been studied. This is be- 
cause immediately after the quench the distribution func- 
tion of the fermions is highly broadened, thus creating an 
enlarged particle-hole continuum. As a consequence for 
attractive interactions either no solutions are found, or 
only damped solutions that lie in the particle-hole con- 
tinuum are found. Further, if the current in the athermal 
steady state is larger than a critical value, then even the 
undamped mode for repulsive interactions vanishes in the 
long- wavelength limit. 

These results are consistent with the ones obtained in 
Ref Hi] where it was found that as a result of a quan- 
tum quench and mode-coupling, a Luttinger liquid de- 
scription is replaced by a low energy effective theory of 
thermal bosons with a finite lifetime. Indeed, in equi- 
librium, the undamped collective modes obtained from 
RPA can also be described as a Luttinger liquid The 
RPA analysis shows that for attractive interactions the 
collective modes lie in the particle-hole continuum and 
are therefore overdamped. The analysis of the present 
paper thus allows one to interpret the generation of the 
friction that was put in evidence by the RG analysis of 
Rcf 21 as due to a generalization to an out of equilibrium 
case of Landau damping. 



The paper is organized as follows. In section [IT] the 
models that will be studied and the notations and con- 
ventions are defined. In section IIIII the RPA analy- 
sis where the quench is from the gapless phase of the 
transverse-field Ising model to the XX model is consid- 
ered. In section HVl RPA involving fermions in a current 
carrying steady state is presented, and in section |V] we 
summarize our results. 



II. MODEL 

Below we describe the two different quenches which 
lead to non-equilibrium steady states without (sub- 
section [HX]) and with f sub-section III Bp currents. 



A. Quench from ground state 

The XY spin chain in a magnetic field is defined as 

j 



(1) 



where 7=1 corresponds to the transverse-field Ising 
model. The XY model has been extensively studied)^"— 
and its equilibrium properties are well understood. It is 
also a popular modeU^i^i'^^ for studying non-equilibrium 
situations due to its simple mapping to free fermions. 
Writing this Hamiltonian in terms of Jordan- Wigner 
fermions 



Sj' — c'j exp 



Sj = exp 



n<j 



n<j 



= c]c 



3-3 2- 



(2) 

(3) 
(4) 



we obtain 



(5) 



This is diagonalized by a Bogoliubov rotation^^ 

H^^J2'^vUk, (6) 



where 



= — Jsign ( cos ^ ^ 'Jj \ (cos ^ 'j] +7^ sin^ fc(7) 
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and 



with 



Cfc 



V-k 



cos 6k = 



sm tik 



\cosk-{h/J)\ 



(cosfc - -y) + 7^ sin^ fc 
sign(cos k — {h/J))j sin fc 

(cos fc - j)^ + 7^ sin^ fc 



(8) 

(9) 
(10) 



Here Cj — -^^k^^^'' ^k- The ground state is obtained 
by occupying all modes with negative energy. We will 
be interested in the special case of7 = /i/J=l, where 
the system is critical. ■^•^ The ground state is defined by 
Vk |$o) — for all fc, as — 2j|sin|| is always non- 
negative. 

Given this initial state, we perform the quench by sud- 
denly switching off the anisotropy 7 and magnetic field 
h. The subsequent time evolution is due to the XX- 
Hamiltonian, 



H 



XX 



qV qv 



E 



(11) 

(12) 



where efe = —J cost and the Ck are the momentum- 
space Jordan- Wigner fermions defined above. At long 
times after the quench, the system approaches a diagonal 
ensemble. To see this, note that immediately after the 
quench, the following quantities are fixed by the initial 
state 
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{ckC-k)o 



,2 ffc 



vlvk 



- sin' ^Lkvlk) (13) 



. . Cfc Uk 

I sm — cos — 
2 2 



rikTll 



vlkV-k 



(.14) 



Time evolution of the c-operators takes the simple form 



Cfe(t) 



Cfc. One finds 



4(0cfc(t) 



4 Cfc 



(cfc(i)c_fc(t)) = e-2"'=*(cfcC_fe)o, 
When averaged over long times, 
1 



tb — ta 



dt {Ckit)c_k{t)) 



(15) 

(16) 
(17) 



for {tb — ta) 00, due to the rapidly oscillating exponen- 
tial. Thus we obtain a diagonal ensemble in the long-time 
limit with a highly broadened momentum distribution 
given by 



4 Cfc 





. k 




sm — 


) 
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(18) 



Note that, as we discuss at the end of this section, such 
an approximation is not necessary and one can retain 
the oscillating modes. It however considerably simplifies 
the expressions to explicitly eliminate them. By equat- 
ing the above non-equilibrium distribution function to a 
Fermi function, one may define a (momentum dependent) 
effective temperature^ 



^cff — 



-Jcosk In 



1 



sm ; 



1 - I sin 11 



(19) 



where one notes T^^ < for fc < 7r/2. As we shall see 
later, since the system is out of equilibrium, this temper- 
ature is not universal, but depends on the quantity being 
studied. 

Once the above steady state has been reached, we con- 
sider the effect of nearest-neighbor Ising interactions in 
the XX Z model 



Hf^H 



XX 



(20) 



where we assume that was switched on very slowly, so 
that in the absence of a quench, the fermions will evolve 
into the ground state of the XX Z chain. The effects of 
this interaction term will be treated within the RPA. 
The basic fermionic Green's functions defined by 



Gfik;t,t') 
Gf{k-t,t') 



-^e{t-t'){{ck{t),cl{i;)}), (21) 



are found to be 



Gfik;t,t') 
Gf{k;t,t') 



Ck{t),cl{t'] 



-ie{t - ty^""-^'-*") 
k 



(22) 



-iefc(t-t') 



(23) 

sign(e),). (24) 



Within the RPA, the particle-hole bubbles ari 



n«(l,2) = ^[Gf(l,2)Gf(2,l) 

<K 



,35 



Gf(l,2)G^(2,l)] 



(25) 



n^(l,2) = ^[Gf(l,2)Gf(2,l) + Gf(l,2)G^(2,l) 

+ G;(1,2)G«(2,1)]. (26) 
which in frequency-momentum space are given by, 
dk dfl 



271 2n 



[Gf{k + q,uj + n)Gf{k,n) 



+ Gf{k + q,Lu + n)Gf{k,n)] (27) 

+ Gf{k + q,uj + n)Gf{k,n) 

+ Gf{k + q,uj + n)Gf{k,n)]. (28) 
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The collective mode dispersion is defined by the roots of 
the complex dielectric function;^ 



eRPA((7,a;,) = 1 - VqIl^{q,LL>q) = 0, 



(29) 



where we neglect the g-dependence of Vq and take it to be 
Vq = J'^ = Vq. The RPA analysis is given in Section Hill 
As mentioned above, it is not critical to work with 
the diagonal ensemble. If we had retained the full time- 
dependence in Eq. (|16p , the integration over the internal 
k variable in the evaluation of the RPA bubbles would 
result in terms that decay with time. Since we are ulti- 
mately interested in the long-time limit, these contribu- 
tions are not important for us. 



B. Quench resulting in a current-carrying state 
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FIG. 1: Initial c-fermion distributions for (a) no current (fco = 
0), (b) non-zero current (fco = 7i"/2), and (c) majcimum current 
(ko = tt). Note the sharp discontinuity for the case of non- 
zero current. 



We will also be interested in how the collective dy- 
namics change when the athermal steady state is char- 
acterized by a net current. This is generated by adding 
a Dzyaloshinskii-Moriya interaction term"^- to the XY 
model. Hi ^ Hi + Ho m , where 



Hdm - -a XI [^P'Ui " ^j^j+i 



(30) 



This new Hamiltonian is diagonalized by the same Bo- 
goliubov rotation"^^'"^^ that diagonalizes the pure XY- 
model. In the isotropic case {XX chain), this can be 
interpreted as a spatially dependent, physical rotation of 
the spinSi ^^'^" For the more general anisotropic chain, the 
spectrum is similarly modifie d^^'^^ 



H,{X) = 



k 



ft 



A sin k. 



(31) 
(32) 



with given in Eq. ([7]). A has the effect of raising the 
energies of states with fc < while lowering the energies 
of modes with positive k. The occupation number is 
now nonzero for A > J, and the 77-fermion occupation 

e{k)9{ko - k), with 



IS Uk 



„ -^1 J 
ko — 2 cos — . 

A 



(33) 



We will see in section IIVI that the presence of this 
non-zero "Fermi momentum" will give rise to multiple 
damped modes within the particle-hole continuum that 
are not present for the zero-current steady-state. Fur- 
thermore, above a certain critical filling factor the single 
undamped collective mode will cease to exist. 

The asymmetry in momentum space drives a current 
in the modified ground state given by 



(j„) = JIm^5++i5„ 
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(A< J) 
(A> J) 



(34) 



It should be noted that this operator can be interpreted 
as the current operator only within the XX-model where 
the total magnetization commutes with the Hamiltonian. 

Performing a quench where A, h and 7 are switched off 
allows this state to evolve under the XX model, obtain- 
ing a non-equilibrium momentum distribution 



clck 







. k 




sm — 






2 



e{k)e{h) - k) 



e{-k)e{ko + k) 



1 



sm ■ 



(35) 



and a current given by Eq. (j34[) . The distribution func- 
tion for several different current strengths is shown in 
Fig.ffl 

After the decay of transients (cfeC_fe), (^cj^clj,^, we in- 
vestigate the collective modes by employing the RPA 
analysis outlined in sub-section III Al The RPA requires 
knowledge of the single-particle Green's functions. The 
presence of a current does not affect the retarded Green's 
hmction (Eq. ((23)) ). but modifies the Keldysh Green's 
function as follows 



Gf{k-t,t') = -ze-"'^^*-*') 
- (nfc - n_fe)] , 



k 



(1 - Uk - n^k) 

(36) 



where = 6{k)9{kQ — k) is the occupation number of 
the ry-fermions in the initial state. 



III. RPA FOR QUENCH FROM GROUND 
STATE 

In this section, we investigate the effect of interactions 
on the athermal steady state (sub-section III Ap obtained 
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from quenching from the ground state of the transverse- 
field Ising model. The RPA particle-hole bubbles are 



dk 
2^ 



cos6lfe(l - 2nk) 



{u} + iS) + Cfe - ek+q 
cosOk+qjl ~ 2nk+q) ' 
(a; -I- iS) + ek- ek+q_ 



(37) 



of a; > 2 J sin | there are no roots to the argument of the 
delta function in the Keldysh component, and we have 
{q, uj) — 0. With n^(g, a;) purely real, this regime lies 
outside the particle- hole continuum. In sub-section IIII CI 
we will demonstrate the existence of an undamped col- 
lective mode lying just above the particle-hole continuum 
for repulsive interactions only. 



n^(g,a.) = ^(27r) j —5{i, + ek~ek+q) 
x[cos0fccos0fe+,(l - 2nfe)(l - 2nk+q) - 1]. (38) 

For the case of interest, 7 = /i/J = 1, we have cos6'fc — 
I sin 1 1, and the distribution of the 77 fermions Uk — 

There are some basic symmetries of the po- 
larization bubbles that are worth mention- 
ing. Firstly n^>^(q,w) = n^'^(-q,w), while 
Re[n^](g,-w) = Re[n^](g,w),n^((7,-a;) = n^((7,w) 
and Im[n^(q, — w)] = — Im[n^](g, w). Therefore in what 
follows we will assume g > 0,a; > 0, and the results for 
the other regimes can be extrapolated from the above 
symmetries. 

There are two regimes which we will study separately. 
One is w > 2 J sin | where Imp-'^] = 11^ = 0, and the 
other is w < 2 J sin f where Im[n^] 7^ 0, 7^ 0. 



Evaluation for to > 2 J sin | 



In this regime, the integrand contains no poles, and 
the result is purely real: Re [n^(g,a;)] = n^(g,a;) and 
Im [n-'^(q, w)] — 0. One may safely take 5 — >■ to find 



COS ■ 



(2Jsinf)2 





1 + 


sin f 




1 + 


sin j 




In 




— z_ In 






1 - 


sin f 


1 - 


sin j 










2_ 





sm ■ 



2-Ki^uj'^ - (2Jsin§)2 
'1 + ^ 



In 



1 - 



2+ 



z_ In 



(39) 



with 




Vw < 2 Jsin ■ 



1 Vw > 2 Jsin I (40) 



Note that our convention is to place the branch-cut of the 
logarithm on the negative real axis. In this same regime 



B. Evaluation for cj < 2 Jsin | 

For w < 2Jsin| the integrand generically contains 
poles. We extract the real and imaginary parts in the 
usual way by writing, 



dk 



f{k) 



271 UJ + iS + Cfe - Ek+q 
dk 



dk. 

2^ 



V 



fik) 



2tt 



f{k)5{uj + ek - ek+q) 



Re [n^] -I- ilm [IL^ 



CO + Ek — Ek+q 

(41) 
(42) 



where V denotes taking the principal value of the inte- 
gral. We obtain 



Rc [n^iq,Lo)] = 



cos ■ 



2^^(2Jsin§)2-u;2 





1 + 


sin f 




1 + 


sin j 








— z^ In 




j 


1 - 


sin 1 


1 - 


sin j 










Z- 





27r^(2Jsinf)2 
1 



In 



1 



In 



(43) 



with z± defined in Eq. (^(1)) . The results for the imagi- 
nary part, Im [n-'^] = ^^—^p—, and the Keldysh compo- 



nent subdivide the particle-hole continuum into two sub- 
regions, < w < 2Jsin^ I and 2Jsin^ | < w < 2Jsin |. 

For < w < 2 J sin^ | one finds 



Im [n«(g, 



-1 



2^(2 Jsin §)2 



sm ■ 



— sm 77 

2 2Jsin^ 



2^(2 Jsin §) 

1 



cos sm 



UJ 



2 Jsin ■ 



!) 

(44) 



(45) 



whereas in the region 2Jsin^ | < a; < 2Jsin|, the result 
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Im [n^(g, 



2^(2Jsin 1)2-^2 



sm — sm 
4 



2'"' 2 J sin ^ 



2'''' 2Jsina 



2^(2Jsin 1)2-^2 



COS 2 

2 



(46) 



(47) 



The particle- hole continuum, which is the region in g 
space where Imll^ 7^ is indicated as the shaded region 
in Fig. [5] and compared with the equilibrium (no quench) 
result (inset). The two different shadings refer to the dis- 
continuity in the functional forms of Ini[n]^, 11^ across 
u) = 2Jsin^ |. The consequences of these results are dis- 
cussed below. 




FIG. 2: Undamped mode (dashed line) above the extended 
particle-hole continuum (shaded region). Inset: equilibrium 
continuum and undamped mode. 



C. Undamped mode for u > 2 J sin | and Vo > 



This can be compared to the undamped mode in the 
equilibrium problemj^i^ 



The undamped mode is obtained from the solution of 

1 - Vbn'"((7,w) = (48) 

with Il^{q^uj) given in Eq. (|39|. For sufficiently small 
Vb, we need to identify the points where diverges. 
This occurs for w = 2Jsin(g/2). Thus in the limit of 

with 



w — 2 Jsin ^, we can write z± ~ \ i: [1 ± . , 



(2 J sin f] 



The dominant contribution as- 



suming that (7 — )• (so that cos | « 1) is given by 



27rn/a;2 - (2Jsin§) 



z+ In 



1 + V2 + i-rT^ 

^ 4 J sm ^ 

1 - V2 + Z . 

^ 4,7 sin t: 



z_ In 



1 + ^2- 



eV2 



4 J sin : 



4 J sin 



(49) 



ttJ' 



(52) 



which exists for both attractive and repulsive interac- 
tions. Thus, for repulsive interactions, the obtained 
sound wave is qualitatively similar to the equilibrium 
case, but with a slightly modified velocity of propaga- 
tion, whereas for attractive interactions, no undamped 
modes exist. 



D. Enhanced particle- hole continuum and effective 
temperature 

The highly broadened initial fermion distribution gives 
one way to define an effective temperature in this non- 
equilibrium state (cf. Eq. (fT9|)). By analogy with the 
equilibrium properties of the particle-hole bubbles, one 
may define an effective temperature in terms of the col- 
lective degrees of freedom. Since the system is out of 
equilibrium, this temperature will in general depend on 
uj,q and also the chosen correlation function 



Due to the branch-cut in the logarithm (chosen to be on 
the negative real axis), in the limit e — )■ 0, the above 
expression can be further simplified to give 



U\q,u;) 



2Juj^ - (2Jsin§)2' 



(50) 



Thus from Eq. (j48|) we find a single undamped mode 
provided Vq is positive with a dispersion 



UJg^eiVo)J\q\\ 1 + 



^0 



32J2- 



(51) 



n^{q,io)-n^{q,io) 



= tanh 



n^(g,^) ^^""V2T,'ff(g,w) 
For small frequencies, this ratio yields 



2 J sin f 



(sm- 



(53) 



cosjj, (54) 



so that for — > 0, we obtain an effective-temperature 

J sin I 



(55) 
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We argue below that T^ff{q,u! — 0) is responsible for 
smearing out the particle-hole continuum in much the 
same way that temperature in an equilibrium system 
does. To see this recall that the particle-hole continuum 
represents the region in the {q, w)-plane where a collec- 
tive mode of frequency uj and wave-vector q is unstable 
to decay into single particle-hole excitations. The upper 
and lower continuum boundaries in equilibrium at zero 
temperature are given hy^^^ 



^L{q) = J sing, 
Louiq) = 2 J sin |. 



(56) 
(57) 



A simple argument makes these boundaries plausible: 
consider the energy of a single particle-hole excitation 
which is created by removing a particle of momentum k 
and creating a particle of momentum k + q, 



uj{k, q) = ek+q - Efc = 2 J sin | sin 



(58) 



This excitation energy depends not only on the momen- 
tum of the excitation, but also on the momentum of 
the original particle, k. For a half-filled band at zero 

temperature, the occupation of fermions is ' 



C^Ck 

6* (-1 — I A; I). The only momenta available for hole cre- 
ation are those with |fc| < kp = Because the cosine 
dispersion has maximal slope at fc = ^, the maximum 
excitation energy occurs for a given q with ^ = f- — f ■ 
The smallest excitation energy for a given q at zero tem- 



perature occurs at k 



or k 



Thus 



(i 



2 J sin ■ 



1(9) = ( 2'*? ) = ^ ^^^1- 



(59) 
(60) 



which are just the upper and lower boundaries of the 
particle-hole continuum (inset Fig. Now, consider 
lowering kp- Excitations of smaller energy for a given 
q are now possible, and in the limit fc^ — 0, we have 
'^min — ^ 2Jsin^ |. The result is the same if one considers 
the opposite limit oi kp ^ n at zero temperature. 

In the present non-equilibrium situation, we find a 
particle-hole continuum (Im [n^(g, oj)] 7^ 0) that ex- 
tends below this lower-bound all the way to a; = 0. 
A finite temperature is known to smear out this lower 
boundary,^** due to the smoothing out of the zero- 
temperature step function for the occupation proba- 
bility. It is interesting to note that the expressions 
for Im [n-'^(g,a;)] andn^(g,a;) are actually continuous 
across the line uj = 2 Jsin^ |, with discontinuities appear- 
ing in their derivatives. In Fig. [5] we plot the undamped 
collective mode dispersion with the particle-hole contin- 
uum represented by the shaded region. The two different 
shadings are separated by the line u) = 2Jsin^ |. The 
analogous plot for the equilibrium situation is shown in 
the inset. 



IV. RPA FOR CURRENT CARRYING STATE 

We now apply RPA to study the current carrying non- 
equilibrium steady state described in section III Bl The 
Keldysh component of the fermion Green's function is. 



(61) 



iGf{k;t,t') = [cos6'fe(l-nfc-n_fc) 
- (nfe-n_fe)]e-"^(*-*'). 



while the retarded Green's function is given in Eq. 
and Uk = 9{k)9{ko — k). Eq. ([6T|) implies that the dis- 
tribution function for the Jordan- Wigner fermions in the 
current carrying post-quench state is not only broad as 
for the zero current case, but is also asymmetric in fc, 
with sharp discontinuities superimposed on it (see Fig.[T]). 
Thus we will find that as for the zero-current case, the 
particle-hole continuum here too is broadened (extend- 
ing everywhere below the line a;,nax = 2 J sin |), while 
the sharp structure in the distribution gives rise to some 
discontinuities in the expression for Im[n^], and the ap- 
pearance of additional damped modes. 

The particle-hole bubbles are now given by 



n^(9,^) 



1 r dk 
'2 J 2^ 

(jik 



cos 9k 



^fc(l - nk - n_fc) 
LLi + iS-2Jf 

- n^k) 



sin I sin {k + §) 



UJ + i5 - 



2Jsin|sin(fc + I 

C0s9k+q{l - Uk+q - n^k-q 

i(5-2,/sinf 







UJ + iS - 



I sin (fc + f ) 

n-k-q) 

2Jsinf sin(fc + f) 



(62) 



Ti^{q,uj) = -{2n) 



dk 
2^ 



5{uj + ek- ek+q) 



X {[cos 6*^(1 - Uk - n^k) - {nk - n^k)] 
X [cos6'fc+g(l - Uk+q - n^k-q) - {nk+q - n^k-q)] 

-1}.(63) 

As before two regions appear, one where w > 2 J sin | 
for which Imfll^] = 11-^ = 0, and the second being w < 
2 J sin I where a particle-hole continuum is found to exist. 
We discuss these two regions separately. 



A. Evaluation for u> > 2 J sin | 

In this regime, as before, the result is entirely real 
and we let (5 — > 0. We find it convenient to write 
Y[R ^ jj(i) _|_ jj(2)^ where n^^' depends on the Bogoli- 
ubov angle, cos 9k ~ |sin-||, while 11^^^ contains the 
rest. As before it is convenient to summarize the symme- 
tries of the polarization bubbles. We find n(^^(— g, w) = 
n'^^'(g,a;), however due to current flow, n(^^(— (7,cj) — 
-n(2)(q,tj). Similarly, Re[n(i)](g, -w) = Re[n(i)](g,a;), 
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while Re[n(2)](5,-w) = -Re[n<-'^^]{q,uj). In the discus- 
sion that follows, we take g > 0,^ > 0. 

We find, 

nW(g,c.) = - 



T sin 



2 



47ri 



T cos 



(2Jsinf)2 



2 



(2J sin 1)2 

fco 

2 



sin 



cos 



2 



(P4) 



Re [n«(g,a;)] = Re[n(i)] + Rep^^)], where 
Re[nW](g,c.) = - 



47r^(2 J sin 1)2-^2 



X < J" sin 



2 4 



sm ■ 



47r^(2Jsin§)2 -w2 



X J" cos 



where 



J- sin 



J" cos 



2 4 



fco g 
2 4 



{P8) 



where 



T(z) — z+ In 



— z_ In 



(65) 



J^(z) = z+ In 



and 



1 - ^ 



z_ In 



(69) 



In the limit /cq — ^ 0, we recover the results of section lllll 
The remaining terms can be collected as 



n(2) = 



47r 



r.fco+9/2 i'q/2 
lko-q/2 J-q/2 



X dk 



Ld~2J sin I sin k uj + 2J sin | sin k 
1 



(66) 



47rz^a;2 _ (2Jsinf)2 
1 



In 



1 



t^tan(^ + j) 
4J sin §2^1 

c^tan(^ + j) ' 
4,7 sin "1 2^ 



In 



"tan(»-j) 
' 4Jsin ^z'l 

"tan(^-j) 



In 



2 In 



, ^tan(M + |) 
^ 4Jsinf2^ 

-, ^tan(^+'j)' 
4Jsinf2^ 

u; tan( ^ ) 



4.7 sin f 
a; tan( ^ ) 
4.7 sin ^z'i 



In 



2 In 



4 J sin "1^^ 



4 J sin '. 



cjtan(^-; 
4 J sin f 

_ '^tan( j) 
4 J sin f zi 
tan(^) 
" 4Jsin ^z'-' 



H67) 



The consequence of the above expressions for n(^'2) 
will be discussed in section ITV CI 



In 



In 



Re[n(2)](g,w) = 

„tan(^l+|) 



47r^(2Jsin§)2-a;2 



1 + 



4 J sin 1 



_ a;tan(^ + j) 
4 J sin f z^ 

^tan(M + j) 
^ 4.7 sin 2 2^ 



^tan(^ + f) 
r TTTTT— oTS 



- 2 In 



1 + 



4 J sin ^z'i 



4 J sin : 



oj tan 



(I) 



4 J sin -^z- 



-In 



In 



2 In 



1 + 


w tan(^-^ — -^^ 


4 J sin 


1 - 


"tan(^-l) 


4 J sin 2z^ 


1 + 


<,tan(M_|) 


4 J sin f 2^ 




1 - 


^tan(^-l) 


4Jsin f 2^ 




1 + 


w tan( ^ ) 




4J sin 


1 


1 - 


u tan( ^ ) 




4Jsin f z^ 





(70) 



We do not give expressions for Im [n^] , 11^ as the 
boundaries of the particle-hole continuum are the same 
as in Fig. [2l though there are additional discontinu- 
ities within the continuum besides the one along w — 
2Jsin^ |. Instead in the subsequent sections, by study- 
ing the divergences in Re [11^] we will identify a single 
undamped mode for repulsive interactions and fco < 7r/2, 
and several damped modes for both repulsive and attrac- 
tive interactions. 



B. 



Evaluation for uj < 2 J sin | 



In this regime, the real part of X\^\q^ijj) is given by 
the principal value of the integral in Eq. (15^ . Writing 



In this section we demonstrate that an undamped col- 
lective mode survives provided the current is not too 

large. As before, define e = \J lo"^ — (2Jsin |)^. In the 
limit g — > and for small fcoj the most divergent terms 
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in n^((7, uj) are 

n«(g,w) ^ 

+ T 



sm ■ 



4^J(2Jsinf) 



T 



ko 



cos 



fco 
2 
sin ■ 



(71) 



27riV2^w2 - (2Jsin|)2 



fco ^ .e\/2cos^ 



1 - \/2 cos — + i . 

2 2 J sm ^ 



In 



1 - V2i 



fco 



.e%/2cos^ 
^ 2 J sin § 



(72) 



The above expression shows that provided fco < -I , which 
corresponds to the logarithms having a branch-cut, we 
obtain 



n«(q,c.) 



2^a;2 - (2J sin 1)2' 



(73) 



Thus by setting 1 — VqII-'^ = 0, we recover the same dis- 
persion as in the absence of current 



M^^ + ^0{Vo)d{7T/2-ko) (74) 

Thus, the undamped mode is unchanged for a current 
which is below the threshold value of fco < tt/2. On the 
other hand, for currents larger than this value (fco > 7r/2) 
and for q ^ kg, no undamped modes exist. 



In this regime, all modes are damped. We iden- 
tify these damped modes by looking for solutions to 
1 - VoRe[U^]{q,ujq) = 0. For Vq 0, all we need to 
do is identify where Re [^^] ±oo. Then positive di- 
vergences correspond to damped modes with repulsive 
interactions, while negative divergences correspond to 
damped modes with attractive interactions. 

Upon examining Eqns. (j68p and (I70|) . we find logarith- 
mic divergences in Re [n^((j', w)] along the characteristic 
lines (for uj,q > 0) 



uji{q) ^ 2 J sin^ |, 

uj2{q) = 2 J sin I sin fco -I- I) 



^iil) = ±2 J sin I sin ( ^ - fco 



(75) 
(76) 
(77) 



Note that W2,3 coincide with the characteristic lines in the 
equilibrium problem with an arbitrary Fermi momentum 
fco4i In the equilibrium problem, these lines represent 




FIG. 3: Characteristic lines along which Re [n^((/, oj)] di- 
verges, giving rise to damped modes for repulsive (solid lines) 
and attractive (dashed lines) interactions for fco = f • Only 
the top line to = 2 J sin | corresponds to an undamped mode 
existing above the particle-hole continuum. 



boundaries across which Im [ll-^] undergoes a jump dis- 
continuity. This is also the case here, though we will 
focus our attention on the behavior of the real part. 

One finds Re [ll^] +oo along the line w^(g) for 
2fco < g < TT and along the line ui2{q) for g < tt — 2fco. 
These correspond to damped collective modes for repul- 
sive interactions. Furthermore, Re [n^] — > — oo along 
the line wi(g) for all q S (0, tt), along uj^{q) for q < 2fco, 
and along u}2{q) for g > tt — 2fco. These negative diver- 
gences represent collective modes created by attractive 
interactions. We plot these characteristic lines in Fig. [3] 
and indicate whether the mode exists for attractive or 
repulsive interactions. 

Such damped modes are usually considered physically 
uninteresting^! compared to any undamped excitations 
in the system, as the damping makes these modes ex- 
perimentally unobservable. The divergences in n^(g,cj) 
that give rise to these damped modes are of a different 
nature than those giving rise to the undamped mode. To 
see this consider the case of w ~ —2J sin | (sin f ± e) 
for small e. The dominant contribution is given by 



Re [n«(g,w)] 



21og^ 
e 



(78) 



where Aq is a g-dependent factor. 

Solving 1 = VoRell^ [(g,a;g)] to leading order in e, one 
finds 



Re \u}„ 



2 J sin 



— sm — ± AnC 
2 V 2 ^ 



(27rJsinq)/|Vo| 



(79) 



where g > is assumed. In general, for any of the char- 
acteristic lines described above, the ansatz w = Wc ± 
2Jesin| leads to a divergence of the form n^(g, cj) ^ 
± log i . Each logarithmic divergence corresponds to two 
damped modes lying exponentially close to each charac- 
teristic line. 
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One may study the modes near lo « 0^2(9) in a way 
similar to our analysis for the modes near uji. For u « 
0J2{q) ± 2Jesin|, one finds 



Re[n«(g,c.)] 



1 + sin^ A' , ^ 

SttJ sm I |cos (fco + I e 



which gives rise to damped modes for repulsive interac- 
tions (Vb > 0) when g < tt — 2fco with 



Re [ujg\ ~ 2 J sin ^ |^sin + /cq 
± A'q exp 



SttJsIu I I cos (fco 



V?, 1 



- sm ■ 



(81) 



The results for the other characteristic line is similar 
and we do not discuss it further. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have applied the RPA to study the ef- 
fect of weak Ising interactions in a non-equilibrium steady 
state of the XX Z spin chain. This non-equilibrium state 
was created in two different ways. One is by quenching 
from the ground state of the transverse-field Ising model 
at critical magnetic field to the XX-model. The sec- 
ond was to modify the Hamiltonian before the quench by 
adding Dzyaloshinskii-Moriya interactions. This had the 
effect of creating a current carrying state. 

The RPA for both the steady-states shows the exis- 
tence of a single, undamped, collective mode for repulsive 
interactions which is qualitatively similar to the sound 
mode in equilibrium, but with quantitative changes to 



the mode velocity (c.f. Eq. ([51]) '). However if the cur- 
rent is larger than a threshold value, this undamped 
mode ceases to exist in the long-wavelength limit (c.f. 
Eq. (TMl) ). The primary effect of the quench is to give rise 
to a highly broadened distribution function (c.f. Fig. [l] 
Eq. ([35]) ) which results in an enhanced particle- hole con- 
tinuum. The boundaries of the particle-hole continuum 
are shown in Fig. [51 Thus for attractive interactions ei- 
ther no modes are found for the first steady-state, or some 
damped collective modes are found for the steady-state 
with current. 

These results, and in particular the generation of a 
finite friction due to an out of equilibrium situation, 
are rather generic and do not depend on the details 
of the non-equilibrium steady-state. Further, the up- 
per boundary of the particle-hole continuum occurs at 
Wmax = 2 J sin I and is related to the fact that the sys- 
tem is on a lattice, and therefore the excitations have a 
maximum velocity. If instead a quadratic dispersion for 
the fermions is adopted, then there is no upper-limit to 
the velocity of excitations. This together with the fact 
that immediately after a quench, the Fermi-distribution 
is very broad with no well defined kp will further en- 
hance the upper boundary of the particlc-holc contin- 
uum, damping even the mode with repulsive interactions. 

An important future direction for research is to ex- 
plore how these results change when an explicit time- 
dependence on is introduced. In particular it is im- 
portant to understand how slowly has to be turned 
on in order to recover the results of this paper. 
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